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ABSTRACT 

We analyze the photometric short-cadence data obtained with the Kepler 
Mission during the first eight months of observations of two solar-type stars 
of spectral types G and F: KIC 11395018 and KIC 11234888 respectively the 
latter having a lower signal-to-noise ratio compared to the former. We estimate 
global parameters of the acoustic (p) modes such as the average large and small 
frequency separations, the frequency of the maximum of the p-mode envelope 
and the average linewidth of the acoustic modes. We were able to identify and 
to measure 22 p-mode frequencies for the first star and 16 for the second one 
even though the signal-to-noise ratios of these stars are rather low. We also 
derive some information about the stellar rotation periods from the analyses of 
the low-frequency parts of the power spectral densities. A model-independent 
estimation of the mean density, mass and radius are obtained using the scaling 
laws. We emphasize the importance of continued observations for the stars with 
low signal-to-noise ratio for an improved characterization of the oscillation modes. 
Our results offer a preview of what will be possible for many stars with the long 
data sets obtained during the remainder of the mission. 

Subject headings: Methods: data analysis - Stars: oscillations - solar-type - individual 
(KIC 11395018, KIC 11234888) - Asteroseismology 
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Introduction 



Helioseis mology has proved to be a powerful tool to directly probe the in terior of 



the Sun (e.g. 



Turck-Chieze 



2001 



Thompson et al. 



2003 



Chaplin & Basu 



20081 ). Thanks 



to years of continuous data, this tool 



improved constraints on so 



Basu et al. 



2009 



Serenelli 



provided a better understanding of the Sun an d 



ar m odels (jChristensen-Dalsgaard et al. 



1996 



Gough et al. 



1996 



20101 ) . However, to establish the broader context and to make 
continued progress on stellar evolution theory, we need to study many stars across the HR 
diagram. Asteroseismology has progressed tremendously over the past decade, driven by 



ground-based observations and several satellite missions (see reviews 



1994 



Bouchy &: Carrier 



2001 



20081 ). The CoRoT missio n (IBaglin et al. 



solar-like oscillations (e.g. 



Bedding et al 



Garcia et al. 



2002 



2009 



Aerts et al. 



2008 



yy 



Brown fe Gilliland 



Bedding fc Kjeldsen 



2006) has provided data on a few stars showi ng 



Mosser et al. 



2009b 



Mathur et al. 



2010b|). 



For the stars observed with a high signal-to-noise r a tio (SNR), the acoustic (p) modes co uld 



be un ambiguously identified (iBenomar et al 



2009 



Barban et al. 



2009 



20101 ) and the signature of magnetic activity could be measured ( 



comp 



2010ah . 



ementing ground-based spectroscopic analyses (IBaliunas et al. 



Deheuvels et 



Garcia et al 



1995 



2010) 



al. 



Metcalfe et al. 



With the launch of the Kepler Mission (IBorucki et al. 



2010 



Koch et al. 



20101 ). we now 



have continuous observations with a longer duration and higher precision, allowing the 
detection of more modes on many more stars. During the first year of Kepler observations, 
five stars were observed continuously for more than 8 months. This is the firs t time 
that such long and continuous observations have been available. 



Chaplin et al. 



(120 IQh 



used one month of data for several stars to demonstrate the asteroseismic potential of 
Kepler. We present an analysis of the 8-month-long time series for two solar-type stars, 
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K 



in 



C 11395018 and KIC 



Campante et al. 



( 1201 lh 



123488qj while the analysis of two additional stars are presented 



According to the values of the stellar parameter s (T e g, logg) given in the Kepler Input 



Catalog (KIC, 



Batalha et al. 



solar-like oscillations ( iChaplin et al. 



2010 



Koch et al. 



20101 ). these stars are expected to exhibit 



20111) • KIC 11395018 and KIC 11234888 have a 
Kepler-band magnitude of 10.8 and 11.9 respectively. Some spectroscopic analyses (Creevey 
et al., in preparation) show that KIC 11395018 is a G-type star with a T e g = 5660 ± 60 K 
(Pinsonneault & An, in preparation). For KIC 11234888, the type has not been determined 
with certainty but it is likely to be a late F star, which has T e g = 6240 ± 60 K. Some 
features in their spectra suggest that they might be subgiants. 

In the following Section we explain how the Kepler data have been processed for 
these two stars. In Section 3, we estimate the global parameters of the stars (background, 
rotation, mean large separation, mean small separation, mean linewidth). We characterize 
the p modes in Section 4, yielding lists of frequencies for both stars. Finally, we discuss and 
compare the results. 



2. Observations and Data Processing 

The data used in this work were collected by the Kepler photometer in the period 
from May 2009 to March 2010, corresponding to the first initial run (Q0), the first roll 
of one month long (Ql) and the next three rolls of three months each (Q2, Q3 and Q4). 
Unfortunately, due to the loss of all the outputs in the third CCD-module on January 

X KIC 11395018 and KIC 11234888 are also known as Boogie and Tigger respectively 
within the Working Group #1 (WG#1) which is responsible for the analysis of solar-like 

stars 
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9, 2010, four CCDs were lost and in particular for these two stars, KIC 11395018 and 
KIC 11234888, we have only 20.87 days of measurements d uring the fourth roll. Thus, 



time series of 252.71 days — with a short cadence of 58.85s (IGilliland et al 



2010)— were 



Kjeldsen et al. 



20 id ) 



available to the Kepler Asteroseismic Science Consortium (KASC; 
through the KASOC databascl on these two targets. In the following, we refer to Q01234 
as the full-length time series. 

After the raw-pixel data were downlinked to the Kepler Science Office, light curves of 
the stars were created by calibrating pixels, estimating and removing sk y background and 



extracting the time series from a photometric aperture (IVan Cleve 



2009|). 



Two types of light curves were available for each star: a raw one suffering from some 
instrumental perturbations, and a corrected one in which housekeeping data have been used 
to minimize those instrumental effects during the Pre-search Dat a Conditioning (PDC) 



Jenkins et al 



20101 ). However, in some cases 



allowing the search for exoplanet transits (e.g. 
part of the low-frequency stellar signal (such as the one produced by starspots) can be 
filtered out. Therefore — inside the WG#1 — we have chosen to work wit h the raw data 



and developed our own method for the corrections (see 



Garcia et al. 



201 lh . 



The light curves were corrected for three types of effects: outliers, jumps, and drifts 
(top panel of Figured]). We have considered as outliers in the datasets the points showing 
a point-to-point deviation in the backward difference function of the light curve greater 
than 3a, where a is defined as the standard deviation of the backward difference of the 
time series. This correction removes ~ 1% of the data points. Jumps are defined as sudden 
changes in the mean value of the light curve due, for example, to attitude adjustments or 
because of a sudden pixel sensitivity drop. Each jump has been validated manually. Finally, 



2 Kepler Asteroseismic Science Operations Center 



http : //kasoc . phys . au . dk/ 
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drifts are small low-frequency perturbations, which are in general due to temperature 
changes (after, for example, a long safe mode event) that last for a few days. These 
corrections are based on the s oftware developed t o correct the high- voltage perturbations in 



the GOLF/SoHO instrument (j Garcia et al. 



20051 ) . We fit a second or third order polynomial 



function to the region where a thermal drift has been observed after comparing several light 
curves from the same roll. Then, the fitted polynomial to the light curve is subtracted and 
we add another polynomial function (first or second order) — used as a reference — which 
has been computed from the observations done before and after the affected region. If the 
correction has to be applied on one border of the time series, only one side of the light curve 
is processed. 

Once these corrections are applied, we build a single time series after equalizing the 
average counting-rate level between the rolls (bottom panel of Figure [T]). A change of 
the average counting rate can happen inside a roll when a change of some instrumental 
parameters occurs. To do this equalization and to convert into parts per million (ppm) 
units, we use a sixth order polynomial fit to each segment. 

As a consequence of these instrumental effects, the light curves from Kepler suffer 
from some discontinuities. For instance, they can be related to the pointing of the 
high-gain antenna towards the Earth -to send all the scientific data every month- or to the 
rolling of the satellit e that needs to maintain a proper illumination over the solar panels 



(IJenkins et al 



2010|). By taking into account these gaps, the Kepler duty cycle in the 



case of the observations of these two stars is limited to 93.45 % of the time. On top of 
that, Kepler encountered some instru mental problems du ring these first eight months of 



measurements that we have corrected (IGarcia et al. 



2011 



). Therefore, the final duty cycle 



of the light curves are 91.36 and 91.34 % for KIC 11395018 and KIC 11234888 respectively. 
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3. Global parameters of the power spectrum density 



The power spectral density (PSD) of both stars is re presented in Figure [21 These have 
been obtained by applying the Lomb-Scargle algorithm (jLomblll976l ) to the "corrected" 
data (see the previous section). In solar-like stars the spectrum at intermediate and low 
frequencies is dominated by the convective movements in their outer layers. 



3.1. Background parameters 



The background of the star is modeled as: 



B(u) = W + au 



8=1 



+ (2irvTi) ai 



where W is white noise, which models the photon shot noise, a and b are two parameters of 
a power law taking into account the effects of slow drifts or modulations due to the stellar 



activity, the instrument, etc., and k is the number of Harv ey laws 



OLi parameterize each of the k Harvey-model contributions (jHarvey 



that are used, <7j, r» and 



1985|). 



In our c ase, equation 



described in 



Mathur et al- 



ls fitte d to the spectrum above 5 /iHz following the procedure 



(l2010al ). We initially assumed only one Harvey profile for the 
granulation. The measured characteristic time r c , amplitude a c , exponent a c , and the 
photon noise level for each of the two stars are listed in Table [TJ 

We have verified that the results are not affected by the filtering of the time series. 
We performed the same fits to non-filtered spectra, where no polynomial function was 
subtracted. We recovered fully compatible results for the granulation background and the 
photon noise level, only the very low-frequency component is affected. 

Finally, we have performed fits by adding an extra Harvey profile to model a faculae 



contribution (IKaroff et al. 



2010). When one ensures that the faculae characteristic time 
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t f is shorter than r c (which is also a free parameter), the fitting procedure converges to 
solutions where Tf vanishes. We recover the previous values for the other parameters. Thus, 
we do not find any traces of faculae in the power spectrum of these two stars. 

The background fits using only one Harvey law are represented in Figure [2] (left panel) 
for KIC 11395018 and KIC 11234888 (top and bottom panel respectively). 

3.2. Rotation period 

To determine the rotation period of the star, we investigate the low-frequency range of 
the PSD. Unfortunately, the standard procedure applied to process the data, as described 
in Section 2, filters out the lower part of the PSD below 1 /zHz. We have computed new 
merged time series where we applied a different high-pass filter by subtracting a triangular 
smoothed lightcurve with window sizes between 12 and 20 days. The influence of the 
instrumental drift in the very low frequency domain can thus be investigated. 

The data from Q0, Ql, and Q2.1 (the first month of Q2) show some instabilities due to 
instrumental effects. For KIC 11395018, we have computed the PSD of the full-length time 
series as well as of two subsets of 126 days each. A zoom of the low-frequency region of the 
PSD below 2 /zHz is shown in Figure [3j In spite of the instabilities present in the first half 
of the time series, for the full-length and the half-length subseries the highest peak is at 
~ 0.32 /iHz (Figure [3]), corresponding to a period of ~36 days j^j^ days, where we use the 
width of the resolution bin of the PSD of 0.0458 /iHz to compute the error bars. A closer 
inspection of the PSD of both subseries also implies that the main power is concentrated 
around this frequency but with a relatively more important contribution of the second 
harmonic ~ 0.7 /iHz. 

The case of KIC 11234888 is different as the Q0, Ql, and Q2.1 data are more unstable 
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than for KIC 11395018 and the PSD of the first half of data is dominated by noise. 
Therefore, we have computed the PSD for time series from Q2.2 (the second month of Q2) 
to Q4, which is shown in Figure HI Here, we see one main peak at 0.42 and a group of 
peaks around 0.60 yuHz. This pattern would correspond to a rotation period of the stellar 
surface between 19 jzf'gg days and 27 jl|; 6 g days. The resolution is of 0.065 /zHz because of 
the shorter time-series analyzed for this star. The fact that we observe several high SNR 



peaks may suggest the existence of differential rotati on on the ste 



ar surface as observed in 



200981 ). We notice that 



similar types of stars showing solar-like oscillations ( IMosser et al.l 
the highest peak is ~2.5 times lower than the one in KIC 11395018. This could be due 
either to a low value of the inclination angle or to a less important surface activity leading 
to a smaller impact of the starspots in the light curves. However, for this star, we would 
need more data to confirm this rotation period. 

Finally, we can say that these peaks detected at low frequency are likely to be of stellar 
origin as both stars were observed by the same CCD under the same conditions and show 
different periodicities. 



3.3. p-mode global parameters 



(Mos ser fc Appourchaux 



20081), QML ([Roxburgh 



2009), OCT fHekke r et al 



20091 ). SYD flHuber et al. 



Several pipelines (AAU (Campante et al.ll2010l). A2Z flMathur et al.l l2010ah. COR 



201 Oh . ORK (tested in 



Bonanno et al 



20091 1) analyzed the eight months of data 



to retrieve the global parameters: the mean large frequency separation ((Az/)), the position 
of the maximum amplitude (f max ), the mean small frequency separation, (602), and the 
mean linewidth of the modes. 



Au is the spacing between the frequencies of modes with the same degree (I) and 
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consecutive radial orders (n) and depends directly on the sound travel time across the star. 
It is a very valuable paramet er as it allow s us to estimate the acoustic radius of the star 



and the mean stellar density (jUlrich 



19861 ). Actually, the large separation is not a constant 
and varies with frequency, as shown in Section 3.3.2 so we also calculate the mean value of 
this quantity over a frequency range, (Az/). 

The mean small separation, (£02), is the mean value of the separation between two 
modes of consecutive radial orders and of degrees I = and 2. This quantity is sensitive to 
the structure of the stellar core and provides information about the age of the star. 



3.3.1. Estimation of the mean large separation and z/ max 

The different pipelines use different methods to compute (Az/) and z/ max : either the 
autocorrelation of the power spectrum or the autocorrelation of the time series. Having 
obtained comparable values for the two stars, we put one set of results in Table [2J For 
KIC 11395018, (Az/) is obtained in the range [570, 1140] /iHz while for KIC 11234888, the 
mean value is calculated in the range [465, 935] /xHz. 

KIC 11234888 has a lower (Az/) and z/ max than KIC 11395018. 

The acquisition of longer time series allows better constraints on these global 
parameters with an improved precision. For KIC 11234888 the addition of Q34 (i.e. Q3 and 
Q4) to Q012 enabled us to make an estimate of (Az/) whereas with Q012 alone, we could 
not estimate it with certainty. 

The (A z /) an d z/ max values obtained are consistent with the relationship derived by 



Stello et al. 



(j2009|), viz. Az/ oc v\ 



,0.77 

max - 



13 



3.3.2. Variation of Av with frequency 



As initially proposed by iRoxburghl ( 120091 ). the autocorrelation of the time series 
computed as the power spectrum of the power spectrum windowed with a narrow filter can 
provide the variation of the large separation with frequency, Au(u). The implementation 
of the method, the definition of the envelope autocorrelation function (EACF) and 
its use as an automated pipel in e for th e spec trum analysis have been addressed by 



Mosser fc Appourchauxl (120091 ) . iMosserl (120101 ) showed that, with a dedicated comb filter 



for analyzing the power spectrum, it is possible to obtain independently the values of the 
large separation for the odd and even degree of the ridges (respectively Ai/ dd and Az/ even ). 
Then, proxies of the eigenfrequencies can be derived by taking the highest peak in the 
region where the mode is expected from the variation of the large separation. 

The values of Az/ even and Au ^ for KIC 11395018 are given in Figure [51 The high SNR 
of the data for this star allows us to use a narrow filter for a detailed EACF analysis. The 
method shows the regularity of the even value Az/ even . It clearly emphasizes the lower values 



and the rapid variation of the large separation Au odd compare d to A*/, 



may expect an irregular echelle spectrum due to mixed modes (IMosser 



2010). 



rom this, we 



Results for KIC 11234888 are given in Figure [6j Compared to KIC 11395018, the lower 
SNR for this star requires the use of a broader comb-filter. However, the method shows 
unambiguously the differences between the even and odd values of the large separation and 
exhibits the low values of Az/ Q dd- Again, the presence of mixed modes is suspected. 



3.3.3. Mean small separation and linewidth 

To obtain the mean small separation and mean linewidth, first the spectrum is 
smoothed over fifteen frequency bins. The spectrum is then folded using the value of the 



-14- 



large spacing determined from fitting the modes so that the curvature of the echelle diagram 
is taken into account so that we could obtain vertical ridges. For each degree I, we stack 
the peaks over seven large separations. This leads to an average or collapsed peak for each 
p-mode degree I. Two Lorentzian profiles were simultaneously fitted to the collapsed peaks 
of the I = and I = 2 p modes. 

The distance between the central position obtained for each Lorentzian with its error 
is taken as the average small separation, (802), while the mean linewidth of the modes is 
computed from the fit to the average I = modes. The non-linear least-square fitting 
procedure gives values for the linewidth and the central position of the collapsed peaks. 
Their associated error bars are calculated with the variance matrix. The 1-a error bar of 
the small separation is obtained by using a normal propagation of errors. 

For KIC 11395018, we obtain a mean small separation of 4.12 ± 0.03 /iHz in the range 
[650, 1000] /iHz and for KIC 11234888, 2.38 ± 0.19 /iHz in the range [550, 900] /iHz. For 
the mean linewidth of the modes, (r), we obtain: 0.84 ± 0.02 /iHz for KIC 11395018 and 
0.86 ± 0.06 /iHz for KIC 11234888. 



4. Characterization of the p modes 



4.1. Fitting methods 



Eleven teams fitted the modes of these two stars. We briefly describe here the methods 
used by the different teams. Further details may be found in Table [3J 

The methods used fall into three basic categories. First, a majority of team s adopted 



Appourchaux et al. 



19981 ). i.e., the 



a maximum likelihood estimation approach (MLE; e.g. 
best-fitting model of the frequency-power spectrum was chosen by maximizing the likelihood 
of that model. The model used was a sum of Lorentzian profiles describing each oscillation 
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mode, plus a background term parametrized as per the description s in Section 3.1 (see, e .g. 



Appourchaux et al. 



Garcia et al. 



2009; 



2008 



Barban et al 



Fletcher et al 



2009 



Benomar et al. 



2009 



Bedding et al. 



2010al : 



20101 ). Because of the SNR of this particular set of data, 



it proved difficult to measure the asymmetry of the modes. Thus none of the teams include 
this parameter in the fit. However further efforts in the future may allow the measurement 



of the asymmetry in this and other Kepler stars. A 
of this technique, the Maximum A Posteriori (MAP; 



ew teams app 



Gaulme et al. 



ied a derived version 



20091 ). where they add 



prior information to the MLE. The most common procedure of fitting is to do it globally 



( I Appourchaux 



20081 ). where the entire frequency range is fitted, meaning the entire set 
of free parameters needed to describe the observed spectrum is optimized simultaneously. 
One fitter, however, chose to fit the modes locally (I = 0, 1, 2 (and 3) together) over one 
l arge separation, sim ilar to what has been done for "Sun-as-a-star" helioseismology data 



( ISalabert et al 



20041 ) 



A second, smaller gr oup of teams employed Bayesian Markoy Ch ain Monte Carlo 



algorithms (MCMC; e.g. 



Benomar 



2008 



Handberg fc Campante 



20111 ). The MCMC 



algorithm maps the probability density function of each free parameter, and as such may be 
regarded as providing more robust estimates of the confidence intervals on the parameters 
than is possible from the standard MLE approach. 

Finally, two teams applied a very different approach that did not involve fitting a model 
to the observed frequency-power spectrum. They adopted eith er a classical pre-whit ening 



of the frequency-power spectrum (i.e. CLEAN algorithm, e.g. 



Bonanno et al 



20081 ). or 



estimation of the frequencies of the highest peaks in the smoothed spectrum. 

Only five teams provided an estimate of the average rotational splitting of the 
non-radial modes, and the angle of inclination of the star (see Section 4.3.1 below). 
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4.2. Methodology to select the frequencies 



A careful comparison of the different fitters' lists of frequencies yielded final lists for 
each star, for use in future modeling work. The pr ocedure we ado p ted to compile the lists 
was a slightly modified version of that described in 



Metcalfe et al. 



(I2010M ) 



At each {n, 1} we compared the N different estimated frequencies and identified a nd 



rejected out lying frequencies based on application of the Peirce criterion (IPeirce 



Gould 



1853 1 



1852 



The Peirce criterion is based on rigorous probability calculation and not on any 
ad-hoc assumption. Hereafter we cite Peirce's explanation of his criterion: "The proposed 
observations should be rejected when the probability of the system of errors obtained by 
retaining them is less than that of the system of errors obtained by their rejection multiplied 
by the probability of making so many, and no more, abnormal observations" The logic calls 
for an iterative assessment of the rejection when one or more datasets are rejected. The 
iteration stops when no improvement is possible. 



Following the work of 
follows: 



Gouldl (118551 ). we have implemented the Peirce criterion as 



1. Compute mean x m and rms a deviation from the sample X{ 



2. Compute rejection factor r from 



3. Reject data if \xi — x m \ > ra 



Gouldl (Il855h assuming one doubtful observation 



4. If n data are rejected then compute new rejection factor r assuming n + 1 doubtful 
observations 



5. Do step 3 to 4 until no data are rejected 
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If the number of accepted frequencies was greater than or equal to the integer value of 
N/2, we included that | n,Z} on a minimal list for the star. This is the main difference from 



the previous procedure ([Metcalfe et al. 



2010bl ). where all the teams had to agree on the 



mode to be put in the minimal list. Inclusion on the maximal list demanded that there be 
at least two accepted frequencies for the mode in question (i.e., we applied a more relaxed 
criterion for acceptance). The additional modes in the maximal list are more uncertain and 
should be taken more cautiously but they can be used for instance to disentangle between 
two stellar models that fit the data. 

With the minimal and maximal lists compiled, the next step involved computing the 
normalized root-mean-square deviation of the frequencies of each of the N teams, with 
respect to the average of the frequencies of modes that appeared in the minimal list, i.e., 
we computed 



^normdev.fc 



i Ki - ('vi)l 2 ,,, 

where k labels the team, v^ kl and a\ x are the frequency and frequency uncertainty returned 
by team k, {y n ,i) is the mean value, over all teams, of the frequency of this mode, and Nk is 
the number of modes fitted by the team k and that belong to the minimal list. The team 
with the smallest cr nor mdev,k provided the frequencies which populated the final minimal and 
maximal lists. 

However, some issues remain with this more robust procedure. For instance, the team 
selected to provide the frequencies of the minimal and maximal lists might not have fitted 
one or several modes that should be in these lists. One way to overcome this issue would 
consist of asking the selected team to reanalyze the data using those additional modes. In 
the future, this procedure will be improved. 
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4.3. Lists of frequencies 

4.3.1. KIC 11395018 

Seven months of data have been analyzed by eleven fitters as described in Section 
4.1. Applying the selection methodology detailed in Section 4.2, we obtain a minimal list 
containing 22 modes, i.e. seven orders, as well as the fitter for whom we obtain the smallest 
normalized rms. 

Eight teams also analyzed the eight-month data sets. Unfortunately, these additional 
21 days of Q4, obtained before CCD 3 was lost, did not allow us to identify any more modes 
in the lists. Table H] gives the minimal and maximal lists obtained for this star. Outside 
the frequency range considered here, some more modes are still visible in the PSD and can 
therefore be fitted. The main issue encountered for these modes was to correctly identify 
their degree before doing the fit. In fact, the fits in these conditions become extremely 
dependent on the guess parameters and the identification given to the code beforehand. 
The existence of possible mixed modes complicates the task and we have decided to take 
a conservative position. These frequency ranges will be explored in the future, when more 
data become available. 



Figure [7| represents the echelle diagram (IGrec et al 



19831 ) folded over 47.9 //Hz with 



the minimal and maximal lists, where frequencies have been obtained with the MAP 
method as a global fitting. We can clearly see the different ridges corresponding to the / 
= 2, 0, and 1. While the ridges for the I = and 2 modes are almost straight, the dipole 
modes present a slope and two modes that have been shifted towards the right side of the 



echelle diagram nea r 763 //Hz and arounc 



"avoided crossing" (lAizenman et al 



968 /iHz. This is a characteristic of the so called 



19771 ). which is related to the presence of a mixed 



mode. This type of mode has the particularity of being supported both by pressure, like 
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the acoustic modes and by gravity, like the gravity modes. They are thus sensitive to both 
the surface and the core of the star. The mixed mode is present simultaneously in p-mode 
and g-mode regions, hence the mixed character resulting in a shift of its frequency. These 
mixed modes are therefore very interesting as their characteristics are heavily dependent 
on how evolved the star is and thus they put additional strong constraints on the stellar 



interior, provic 


ing a \ 


Metcalfe et al. 


2010b) 



2009 



For the fitter with the lowest normalized rms, the mixed mode and the mode around 
667 /xHz have been fitted separately, i.e. with the p modes of the previous large spacing, 
instead of globally. Indeed, if the modes are fitted globally, one assumes that the modes 
have a common linewidth and amplitude at every order ({/ = 0, n}; {I = 1, n}; {I = 2, 
n-l}). Thus the first / = 1 and the mixed modes are treated as "single modes". The fitting 
code fixes all the parameters to the values obtained for the modes previously fitted globally 
and then fits these two modes independently with all the parameters free. 

The maximal list contains three more modes at higher and lower frequency. One 
interesting peak is the one around 1016 /iHz. It is located between the two ridges I = 
and 2, and is identified with a high SNR in the PSD. However, its identification remains 
uncertain as it could be either an / = or an / = 2. Moreover the presence of a mixed mode 
in this region is not excluded and could be responsible for this bumped peak. Note that all 
the modes of the minimal list were also detected by the EACF method of Section 3.3.2. 



Finally, the fitter A2Z RG reported a problem while fitting the sp 
inclination angle, which are parameters that are highly correlated (e.g. 



ittings and the 



Gizon fc Solanki 



2003 



Ballot et al. 



20061 ). Indeed, with the whole dataset, the fit converged to null values, 
while with the Q0123 data, the fit managed to converge to a value for the splitting of 
Us = 0.49 ± 0.14 /iHz and an inclination angle of i = 36.83 ± 10.78°. So far, the reasons 
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for obtaining null values remain unclear but it could either be related to a combination of 
the stochastic excitation, the width, and the splittings of the modes or due to a shift of 
the modes that prevents us from properly distinguishing the individual components of the 
modes. 

On the other hand, the MCMC fitting method used by the fitter AAU gave an 
estimation of the projected rotational splitting v$ = z/ s sin(i) of 0.29 ±0.06 /iHz and an 
inclination angle i > 20°, with a confidence level of 68 % by analyzing Q0123 data. Indeed 



Ballot et al. 



( 120061 ) demonstrated that it is more robust to fit the projected splitting, z/g, 
and the inclination angle instead of fitting the usual combination of splitting and inclination 
angle. The MCMC met hod fits the spectrum including rotational splittings as described by 



Gizon fc Solankil (120031 ) . and using a set of priors. The projected rotational splittings, u t 



s- 



and the inclination angle, i, are free parameters. Uniform priors (equal probability for all 
values) are used for both the inclination and u s . The inclination is searched in the range 
[0-90] degrees while v s is assumed to be in the range [0-2] /zHz. Figure |8] shows the PDF 
(probability density functions) of both fitted parameters. These preliminary results need to 
be confirmed with longer datasets. 



4.3.2. KIC 11234888 

KIC 11234888 has a more complex oscillation spectrum. Furthermore, it has a lower 
SNR, which results in a noisier echelle-diagram (Figure [9]). To compute the minimal and 
maximal lists of frequencies, we had to discard the results of two teams out of the six who 
fitted the modes, as one of them could not identify any modes and the other only fitted 
a very small number of modes, leading to very short frequency lists. As a consequence, 
we have used the results of only four fitters to build the final minimal and maximal lists. 
Applying the method described in Section 4.2, we also selected the fitter with the lowest 
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normalized rms. 

The frequencies obtained by looking for the highest peaks in the power spectrum are 
listed in Table [5j The minimal list contains a total of 16 modes, while the maximal list has 
two more modes. Figure illustrates the minimal and maximal lists overplotted on the 
echelle diagram. 

Like the previous star, KIC 11234888 was analyzed with the Q0123 (which corresponds 
to Q01234 data without the last quarter, Q4) and Q01234 data. By adding 21 days of data, 
two modes are discarded in the lists: an I — 1 at 901 /iHz and an I = 2 at 792 /iHz. 

In Figure [U we can see the ridges for the / = 0, 1, and 2 modes. But for the I = and 
2, the SNR is quite low and the modes are not all obviously distinguishable. For the ridge 
of the I = 1 modes, we notice a very interesting structure showing several avoided crossings, 
due to the presence of several mixed modes. 

We can see that there are some regions in the echelle-diagram where we see some 
power but they are not listed as modes. The power present in the I = 1 ridge was rejected 
according to the Peirce criterion. 

We also notice that the modes selected in the minimal list were also detected by the 
EACF method. 

Unfortunately for this star, no value for the splittings and inclination angle were 
provided. 

We need more data and modeling to better understand this pattern and detect the I = 
and 2 modes with more certainty. 
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5. Discussion and Conclusions 

In this work we analyzed the light curves corrected in a specific way of two solar-like 
stars, a G star, KIC 11395018, and an F star, KIC 11234888, observed by the Kepler 
Mission during ~ 8 months. Unfortunately, both stars were located on the same CCD 
that was on the module that broke in January 2010. Thus, only 8 months of continuous 
measurements were available for asteroseismic investigations. Kepler will continue observing 
these stars for at least 2.5 years. Because these stars will be regularly observed by the 
broken CCD, the data will contain periodic gaps of three months in length. We will thus 
have to deal with these gaps either by using subseries of 270 days, by averaging the power 
spectra of subseries, or by removing the fundamental peak and its harmonics in the PSD 
corresponding to the window function. 

Seven pipelines analyzed these data using different methods to retrieve the 
global parameters of the two stars. We obtained the mean large separation, (Au) = 
47.76 ± 0.99 //Hz and (Au) = 41.74 ± 0.94 //Hz for KIC 11395018 and KIC 11234888 
respectively. As the mean large separation is proportional to the average density of the 
star, we can deduce that KIC 11395018 is more dense. 

KIC 11395018 is expected to be a post-main-sequence star due to the clear avoided 
crossings. Moreover, the presence of larger number of mixed modes in the echelle diagram 
of KIC 11234888 suggests that it is more evolved than KIC 11395018. 

We determined the mean values of the small separation, (S02) — 4.12 ± 0.035 /iHz and 
2.38 ± 0.19 /iHz respectively for KIC 11395018 and KIC 11234888, as well as the mean 
linewidth of the modes, (r) = 0.84 ± 0.02 /iHz and 0.86 ± 0.06 /aHz respectively. Note 
that though these stars have different T eff , ~5660 K for KIC 11395018 and ~6240 K for 
KIC 11234888 , we fi nd similar values of linewidths, which is not expected according to 



Chaplin et al 



(120091 ) who concluded that the mean linewidth of the modes scales like T eff 4 
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An analysis of the convective background for both stars demonstrated that 
KIC 11234888 exhibited the largest granulation time scale of the two stars. No signature of 
faculae was found in either of the stars. 

Several teams (up to eleven) analyzed the time series of the two stars with different 
methods and provided p-mode frequencies. After applying a revised procedure to these 
different sets of frequencies, described in Section 4.2, we selected the modes on which the 
different teams agreed to create the minimal and maximal lists of frequencies for both stars. 
We identified 22 p modes in the range 600 to 1000 /iHz for KIC 11395018 and 16 p modes 
in the range 500 to 900 /iHz for KIC 11234888. These frequencies are in agreement with the 
modes detected by the EACF method. 

The rotational splittings have been tentatively measured for KIC 11395018 by four 
teams. The MCMC gave projected splitting of ~ 0.29 /iHz. Combined with the measurement 
of the rotation period of 36 days obtained from an analysis of the low-frequency region 
of the PSD, we improve the inclination angle constraint to % > 45°. For KIC 11234888, 
which seems to show evidence of a differential rotation, we need more data to constrain the 
splittings and the inclination angle. 

Note that for both stars no I = 3 modes have been detected. Indeed degrees higher than 
2 are difficult to detect because of cancellation effects in spatially unresolved observations. 
Howeve r, clear evidence of I = 3 modes has alre ady been seen in red giants obse rved by 



Kepler (IBedding et al. 



2010b 



Huber et al. 



as in the CoRoT target HD49385 (IDeheuvels et al 



2010h and Co RoT (IMosser et al.ll201ll ) as well 



20101 ). We expect to detect them in 



solar-like stars with longer time series and with higher SNR targets. 

A deeper analysis of the characteristics of the p modes of these two stars, such as 
heights, linewidths, and amplitudes, will be presented in another paper (Handberg et al., in 
preparation). 
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With the global parameters of the p modes and the revised values of T e g, we can 
obtain a first esti mation of the mean densi ty, radius and mass of these stars by using 



the scaling laws ( IKjeldsen &: Bedding Il995f ). which are thus model independent. We 
estimate the mean density of these stars: (p) = 0.173 ± 0.007 g/cm 3 for KIC11395018 and 
(p) = 0.132 ± 0.008 g/cm 3 for KIC 11234888. We remind here the scaling laws: 

R & \(Au)J V3050y V57777 1 j 

M_ f 135 \ 4 /^maxW T cff \ 3/2 

M \(Ap)J V 3050 ) V5777/ ' 1 ' 

For KIC 11395018, we obtain: M = 1.25 ± 0.24 M and R=2.15 ± 0.21 R Q , while for 
KIC 11234888, we have M = 1.33 ± 0.26 M and R=2.4± 0.24 ,R . 

A detailed study based on stellar models, which is out of the scope of this work, is 
now possible. The modeling of this star requires the combination of global oscillation 
parameters, sets of frequencies, and atmospheric parameters (from spectroscopy), in order 
to retrieve very precise values of the stellar parameters, namely mass, radius, and age 
(Creevey et al. in preparation). 

Funding for this Discovery mission is provided by NASAs Science Mission Directorate. 
The authors wish to thank the entire Kepler team, without whom these results would 
not be possible. We also thank all funding councils and agencies that have supported the 
activities of KASC Working Group 1, and the International Space Science Institute (ISSI). 
NCAR is supported by the National Science Foundation. SH also acknowledges financial 
support from the Netherlands Organisation for Scientific Research (NWO). 
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Fig. 1. — Raw flux corrected for the outliers (top panel) and corrected data (bottom panel) 
for KIC 11395018. One point out of twenty has been plotted. The thick red dashed lines 
correspond to the beginning of each roll, while the thin red dotted lines correspond to a 
change of the sub-roll (e.g. Q2.2, Q2.3). 
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Fig. 2.— Left: PSD of KIC 11395018 (top) and KIC 11234888 (bottom) smoothed over a 
boxcar of 5 bins. The dashed line corresponds to the white noise level. The triple dotted- 
dashed line represents the Harvey law. The dotted-dashed line is the contribution of the 
power law. The solid line is the sum of all these contributions. Right: Zoom of the Power 
Spectrum Density of KIC 11395018 (top) and KIC 11234888 (bottom) using Q01234. The 
insets represent the power spectrum of the observation window. 
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Fig. 3. — Zoom on the PSD at low frequency of KIC 11395018 oversampled by a factor of 
5. In black the Fourier spectrum of the full-length series, In green, the PSD of the first half 
of the time series and in blue, the second half of the time series. 
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Fig. A.— Zoom on the PSD at low frequency of KIC 11234888 from Q2.2 to Q4 data and 
oversampled by a factor of 5. 
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Fig. 5. — Az/ Cven and Az/ dd for KIC 11395018. The solid lines (resp. dashed lines) correspond 
to the even (resp. odd) ridge; 1-a error bars are given with the same linestyles. We have 
superimposed the proxies of the eigenvalues derived from the method: diamonds correspond 
to I = 0, triangles to I = 1. The inset at the bottom right gives the size of the comb filter 
used for this star. 
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Fig. 6.— Same as Figure [5] for KIC 11234888. 
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Fig. 7. — Echelle diagram of the eight months of data with the minimal list of frequencies 
(filled symbols) and the maximal list (filled symbols plus open symbols) for KIC 11395018. 
The circles correspond to the / = 0, the triangles to the I = 1, and the squares to the I = 2. 

0.02 1 1 1 1 1 1 1 1 1 




Inclination (degrees) v * (jiHz) 

Fig. 8. — Probability distribution of the splittings and the inclination angle for KIC 11395018 
(bottom left panel). Dark colors refer to higher probability values compared to light col- 
ors. The projected rotation is given in the right panel while the inclination angle is represented 
in the top left panel. In each of these two panels, the red line corresponds to the median 
of the distribution and the blue ones indicate the 68% confidence interval around it. The 
dashed line in the contour-plot indicates the resolution in the power spectrum. 
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Fig. 9. — Echelle diagram of the eight months of data with the minimal and maximal lists 
of frequencies for KIC 11234888. Same legend as Figure [TJ 
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Table 1. Granulation and photon noise parameters of KIC 11395018 and KIC 11234888, 
fitted on the power spectrum density using the filtered Q01234 data set. 





To M 


cr c [ppm] 


Oic 


W [ppmVHz- 1 ] 


a (pprrrVHz-^ 6 ) 


b 


KIC 11395018 
KIC 11234888 


698 ± 33 
869 ± 89 


93.1 ± 0.9 
88.3 ± 1.7 


2.02 ±0.06 
1.68 ±0.11 


3.69 ±0.02 
12.75 ±0.06 


1563.28 ± 1142.16 
4130.78 ± 5129.70 


1.84 ±0.33 
1.67 ±0.56 
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Table 2: Global parameters of KIC 11395018 and KIC 11234888 using Q01234 data. 



Star 


Ai> (//Hz) 


z/ max (//Hz) 


(<Jo2> (/iHz) 


(r) (//Hz) 


P rot (days) 


KIC 11395018 


47.76 ± 0.99 


830 ± 48 


4.12 ± 0.03 


0.84 ± 0.02 


36 


KIC 11234888 


41.74 ± 0.94 


675 ± 42 


2.38 ± 0.19 


0.86 ± 0.06 


19-27 



Table 3. Summary of the fitting methods. 



Fitter ID 


Method 


Splittings 


Angle 




Heights 


Linewidth 


Fit mixed mode 


AAU 


MCMC Global 


Free 


Free 


H(Z= 


=1)/H(i=0)=1.48 


Free for 1=0 


No 










H(/= 


=2)/H(Z=0)=0.5 


Interpolation for other modes 
















with 1=0 mode 




A2Z RG 


MLE Global 


Free 


Free 


H(/= 


=1)/H(Z=0)=1.5 


Same for each order 


Yes 




MAP 


Guess 1 fiRz 


Guess 45° 


H(Z= 


=2)/H(Z=0)=0.5 






A2Z CR 


MLE Global 


Fixed to /tHz 


Fixed to 0° 


H(Z= 


=1)/H(Z=0)=1.5 


Same for each order 


No 




MAP 






HO 


=2)/H(/=0)=0.5 


Guess 1 /iHz 




A2Z DS 


MLE Local 


Free 


Free 




Free 


Same for each order 


No 






Guess 0.8 /iHz 


Guess 45° 










IAS OB 


MCMC Global 


Free 


Free 




Free 


Free 


No 


IAS TA 


MLE Global 


Fixed to /tHz 




HO 


=1)/H(Z=0)=1.5 


Same for each order 


Yes 










H(J= 


=2)/H(Z=0)=0.5 






IAS PG 


MLE Global 


Free 


Free 


HO 


=1)/H(J=0)=1.5 


Free 


No 




MAP 


Guess 1 /iHz 


Guess 45° 


HO 


=2)/H(J=0)=0.5 


Guess 2 /iHz 




OCT 


MLE Global 








Free 


Same for all modes 


No 


ORK 


CLEAN 












No 


QML 


MLE Global 


Free 


Free 




Free 


One linewidth per overtone 


No 


SYD 


Peaks in 








Free 




Yes 



smoothed spectra 
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Table 4: Minimal and maximal lists of frequencies for KIC 11395018 in jjHz obtained with 
eight months of data. 



Order 1 = 1=1 1 = 2 



12 




667.05 ±0.22 b 


631.19 ± 1.36 b 


13 


686.66 ± 0.32 


707.66 ±0.19 


680.88 ±0.45 


14 


732.37 ± 0.18 


763.99±0.18; 740.29±0.17 a 


727.78 ±0.30 


15 


779.54 ± 0.14 


805.74±0.13 


774.92 ±0.16 


16 


827.55 ± 0.15 


851.37 ±0.11 


823.50 ±0.16 


17 


875.40 ± 0.16 


897.50 ±0.15 


871.29 ±0.21 


18 


923.16 ± 0.19 


940.50 ±0.15 


918.10 ±0.28 


19 


971.05 ± 0.28 


997.91 ±0.33 


965.83 ±0.23 


20 






1016.61 ±0.73 b ' 



"Mixed mode 

b Part of the maximal list 

c Uncertain identification, could be an 1=0 
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Table 5: Lists of frequencies for KIC 11234888 in /iHz obtained with eight months of data. 



Order 


/ = 


I =1 


I = 2 


11 




506.72 ±0.20 




12 




563.30 ±0.14 


582.84 ±0.21 a 


13 


627.67 ± 0.19 


594.83 ±0.16 


624.65 ±0.18 


14 


669.35 ± 0.16 






15 


711.63 ± 0.15 


686.34 ±0.17 


708.67 ±0.19 


16 


753.64 ± 0.20 


741.10 ±0.18 


751.84 ±0.21 


17 


794.56 ± 0.20 






18 


836.83 ± 0.21 


815.43 ±0.21 




19 


877.80 ± 0.22 a 


855.67 ±0.21 





"Part of the maximal list 



